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ABSTRACT 

We employ tomographic observations of a small region of plage to study the propagation of waves from the 
solar photosphere to the chromosphere using a Fourier phase-difference analysis. Our results show the expected 
vertical propagation for waves with periods of 3 minutes. Waves with 5-minute periods, i.e., above the acoustic 
cut-off period, are found to propagate only at the periphery of the plage, and only in the direction in which 
the field can be reasonably expected to expand. We conclude that field inclination is critically important in the 
leakage of p-mode oscillations from the photosphere into the chromosphere. 
Subject headings: Sun: oscillations — Sun: magnetic fields — Sun: chromosphere 
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1. INTRODUCTION 

The chromosphere appears as a dynamic, thin, barbed red 
ring around the moon about the ti me of total so lar eclipse. 
Early detailed observations by, e.g.. lSecchil dl 877b . have puz- 
zled solar astronomers ever since. We know now that the 
magnetized chromosphere is dominated by "spicules" at the 
limb and their on-disk counterparts ("mottles" and "fibrils", 
lBeckers| [l968). Over 90% of the non-radiative energy go- 
ing int o the outer atm osphere is deposited in the chromo- 
sphere ( Withbroel ll983l) . The chromosphere requires nearly 
one hundre d times the mass and energ y flux of the corona for 
sustenance ( Withbro e & Noveslll977l) and remains the most 
poorly understood region of the outer atmosphere. The com- 
plex dynamic appearance of the chromosphere, at the inter- 
face of photospheric and coronal plasmas where (on average) 
the magnetic and hydrodynamic forces balance, is incredibly 
difficult to interpret unambiguously and is largely responsible 
for the paucity of our knowledge. 

Recent analyses of observations of the chromosphere 
at high spatial and temporal resolution made at multiple 
heights, guided by state-of-the-art numerical simula- 
tions, have begu n to unlock the mysterie s of the dynamic 
chrom o sphere ([De Pontieu et alJ l2004t lHansteen et all 
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&Jefferiesl 120061: Uefferies et al.1 [ 20061: 
2007at iRouppe van der Voort et all 12007 



2007b). Together, these results build a uni 



fied picture of the formation of (at least one type of) spicule 
resulting from the leakage of global p-mode oscillations 
along continuously evolving magnetic field. 

In the traditional view of the solar atmosphere, 5-minute 
p-mode oscillations should not propagate upward because the 
acoustic cutoff is at periods of 3 minutes. However, it has long 
been known that p-mode oscillations propaga te upward in and 
aroun d magnetic flux concentrations (e.g., iGiovanelli et all 
1 1978b . Several mechanisms have been proposed to explain 
this upward propagation of oscillations that have frequen- 
cies below the canonical values for the cut-off frequency in 
the photosphere, e.g., leakage along inclined field lines (e.g., 
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iBel & Leroyll977l:IZhugzhda & Dzhalilovl 19841) . or radiative 
losses in the photosphere (Robertsl l 19831 Because of a lack 
of high-resolution observations and modeling, this issue has 
not been resolved yet. 

Recently, the advent of higher resolution observations 
and modeling have led to renewed interest in this topic, 
with suggestions th at p-mode leaka ge can lead to f orma- 
tion of s picules (e.g..lSuematsvJfT990llDe Pontieu et aT]l200' 
2007a; Hansteen et al. 2006; Roup pe van der Voort et al 
120071: IHeggland et all 120071) and 5-m inute ocillations in 
coronal loops (IDe Pontieu et al.1 120051) . and that it has 
importa nt consequences for the energetics of the atmo- 
sphere (Jefferi es et all 120061) and the d amping of p-mode 
oscillations (Ide Moortel & Rosnerl 120071) . This has caused 
a renewed focus on understanding the causes for the 
propagation of p-mode oscillations into the atmosphere, 
with numeric al models investigating both proposed mech- 
anisms (e.g.. IDe Pontieu et al.1 120041 : IHeg gland etail 120071 
KhomenkoetalJ |2008). While improvements in numerical 
models are important for a better understanding of this issue 
(especially with respect to the role of wave-mode coupling at 
the 0=1 surface), the advent of space-based observations at 
higher spatial resolution allow for new and crucial observa- 
tional constraints which can guide us towards a resolution of 
this issue. In this Letter, we will provide observational evi- 
dence that p-mode oscillations propagate preferentially along 
inclined field at the periphery of plage. 

2. OBSERVATIONS AND REDUCTION 

We employ data sequences taken wi th the space-borne ob- 
servatory Hinode (Kosu gi et al. 12007b. The Solar Optical 
Telescope (SOT. iTsuneta et all 20081: ISuematsu et~aT1 2008; 



llchimoto et alJ2008blShimizu et al.l2008l) was used to observe 
a small area of decaying plage close to disk center (n = 0.1 1) 
on 2009 January 30 using both the Spectropolarimeter (SP) 
and the Narrowband Filter Imager (NFI). A context image is 
shown in Fig.Q] SP was programmed to make repeating raster 
scans of 4.8" x 61.4", with a resolution of 0.16" and an in- 
tegration time of 1 .6 s per slit. NFI was concurrently used 
to observe 5 positions in the NalDi line, at ±168, ±80, and 
m from the center of the line (see Fig. O. The cadences 
of the SP rasters and NFI line scans are about 55 s and 32 s, 




-150 -140 -130 -120 -110 -100 -90 -80 -70 -SO -50 
AT[arcsec] 



FIG. 1. — NalDi V/I image of the region taken at 08:25:02 UT, several 
hours before the start of the sequence used in this Letter. The approximate 
area shown in Figs.[3]and[4]is indicated by the black box. 
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FIG. 2. — The NalDj line (solid line) together with filter transmission 
curves for a perfect Lyot filter using the NFI design at the 5 positions used in 
this analysis (dashed lines), scaled by the NFI Na I Di prefilter transmission 
(dotted line). 



respectively. 

The data were carefully aligned using Fourier cross- 
correlation techniques. First, the NalDi far and inner wings 
were aligned for each time step separately on the basis of 
their absolute Stokes V signal, using the far blue wing im- 
age as reference. Then, using the intensity image, the core of 
the line was aligned to the far blue wing. The sequence was 
then aligned using the Stokes V signal in the far blue wing. 
The NFI wing data was then resampled using cubic spline 
interpolation to match the observing time of the core images. 
Doppler velocities were calculated by fitting parabolas to each 
NalDi 5-pointline scan. 

For the SP data, a linear minus Gaussian function was fit- 
ted to each spectrum of the Fe I 630.25 nm line. The Doppler 
velocities derived from the repeating SP scans were interpo- 
lated in time to match the NalDi core images. Finally, the 
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FIG. 3. — Phase difference between Doppler shift in Fel and NalDi at 
5.6 mHz. Blue and red respectively indicate positive and negative phase dif- 
ference, i.e., upward and downward propagation. Contours indicate the av- 
erage Stokes V signal in Fel at 2.5% intervals. Arrows in the bottom left 
corner of each panel indicate the shift between Na I D] and Fe I using Na I D t 
as reference. Rows are shifted in north by 3, 0, and —3 NFI pixels of 0.16", 
respectively from top to bottom. Columns are shifted east by 3, 0, and -3 
NFI pixels, respectively from left to right. 



SP data was mapped onto a grid matching the square NFI 
pixels of 0.16". This entails scaling and warping of the SP 
data to correct for slit rotation and other distortions as well as 
the application of a fixed offset between the SP and NFI se- 
quences plus the displacements computed for the NalDj core 
sequence in order to correct for image motion. The fixed off- 
set was computed by optimizing the correlation between the 
SP Stokes V amplitude and the Stokes V signal in the NalDj 
far blue wing. The final sequence contains 211 frames. A 
manual check shows an excellent alignment between the SP 
and NFI sequences. 

We employ a Fourier analysis to derive phase differences 
between the Doppler velocities derived from the SP data and 
those from the NFI data. The cadence of the SP observa- 
tions limits the Nyquist frequency to 9 mHz. A cos 2 apodiza- 
tion window is applied to 10% of the sequence. In order to 
reduce noise, we average the phase difference by summa- 
tion of crosspower weighted by a Gaussian with a FWHM 
of 1 mHz centered on the frequency of interest. Propagation 
along slanted rays is sampled by shifting the Na I Di sequence 
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FIG. 4. — Phase difference between Doppler shift in Fel and NalDi at 
3.3 mHz in the same format as Fig. [J] 



by several pixels in each direction. 

While the core of the NalDj line is expected to sample 
higher layers than the core of the Fe I line, it is not immedi- 
ately obvious which Doppler velocity signal is formed higher. 
In order to deduce empirically which is higher, we examine 
the phase difference of oscillations with 2 i -minute periods. 
Upward propag ation is expecte d in the non-magnetic inter- 
network (e.g. jLites et al.|[l993l) . In the majority these areas, 
we find a positive phase difference (i.e., upward propagation) 
if we assume the NalDi Doppler velocity signal is formed 
lower than the Fe I signal. We thus take the Na I Di Doppler 
velocity signal as the reference. The Fe I line is sampled with 
high spectral resolution, which allows for very accurate deter- 
mination of the line core position. Conversely, the relatively 
wide NFI passband, imperfect suppression of sidelobes, and 
the contribution of the NalDi wings in the determination of 
the Doppler velocity cause this diagnostic to sample deeper 
layers. 

3. RESULTS AND DISCUSSION 

3.1. 3 -minute oscillations 

Figure [3] shows the maps of phase difference at 5.6 mHz, 
corresponding to a period of 3 minutes. The central panel 
shows propagation directly upward in the center of the plage 
region. The general pattern remains the same under northward 



shifts, though upward propagation becomes more pronounced 
over a larger area. There is less propagation in the core of the 
plage when sampling in south-west direction. By introducing 
a sufficiently large shift, patches of downward propagation 
appear. 

While the specifics of these results depend on the magnetic 
configuration of the plage, they are not unexpected. Acoustic 
waves with 3-minute periods can propagate in vertical struc- 
tures in the traditional treatment of wave propagation in the 
solar atmosphere. Many previous studies of oscillations have 
shown that this i ndeed happens in both network and internet- 
work areas (e.g., Rutten 1995; Rutt en et al.l20 04). Oscillation 
power is reduced inside the network compared to the inter- 
network, but enha nced in the ar ea immediately surrounding 
the network (e.g., Krijg er et al.l 1200 lb . No obvious phase- 
difference signal that would correspond to these "power au- 
reoles" is present. 

We find much larg er phase difference inside the plage re- 
gion than out side it.TDe Wiin et ail d2005l) analyzed TRACE 
UV continua {Handv et al] | 19991) . and noticed a peak in their 
phase-difference (fe/,,/) diagram at 6 mHz and at small spatial 
scales reaching 110 degrees, and also a corresponding peak 
in their spatially-averaged phase-difference spectrum of net- 
work areas at the same temporal frequency. While those au- 
thors conclude that these features are likely an artifact in their 
data, the analysis presented here shows a conspicuously sim- 
ilar signature while using wholly different observations. We 
thus conclude that the large phase differ ence measured inside 
the magnetic region here and also by iDe Wiin et al.l (2005) 
is not an artifact of the data or processing, but must be at- 
tributed to solar 3-minute waves in magnetic regions. These 
waves must be compressible, since they would not show up 
as an intensity-modulation the TRACE UV continua other- 
wise. Due to the large phase difference, slow mode magneto- 
acoustic waves are a good candidate. Further study of 3- 
minute waves in regions of strong, vertical field through both 
observations and simulations is warranted. 

3.2. 5-minute oscillations 

Figure|4]shows the maps of phase difference averaged over 
a 1 - mHz range around 3.3 mHz, corresponding to a period 
of 5 minutes. In the central panel, only very little propagation 
is detected. None of the panels show significant propagation 
in the core of the plage region. Introducing a shift invari- 
ably results in phase difference indicative of propagation in 
those areas of the plage region where the field is expected 
to diverge in the direction of the displacement. As an ex- 
ample, the left-upper panel shows enhanced phase difference 
in the north-east corner of the plage region. Conversely, the 
bottom row shows increased phase difference on the southern 
edge. Waves with 5-minute periods propagate predominantly 
in those places where the field is inclined in the direction of 
propagation. 

It is possible that the f3 = 1 surface plays an important role 
in p-mode leakage into the chromosphere. The (3=1 surface 
is expected to intersect the photosphere at the periphery of 
plage. If mode coupling is important in p-mode leakage into 
the chromosphere, one would thus expect it to happen pref- 
erentially at the edges of plage, consistent with the current 
results. 

Only weak indications of downward propagation are de- 
tected, mostly at large shifts, e.g., at (-92.5,7) in the upper- 
right panel. Potentially, signals of downward propagation are 
not detected because they are overwhelmed by signals of up- 
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ward propagating waves. If a weaker signal of propagation in 
the opposite direction is present, the phase difference will be 
underestimated. In the case where both signals are of equal 
amplitude, they can be considered as a standing wave with 
zero phase difference. 

Because of the uncertainty in the sampling heights of the 
two diagnostics used, we should consider that in areas with 
strong magnetic field the heights could be reversed as com- 
pared to in internetwork. However, in such a case, the phase 
difference would be negative, indicative of do wnward propa- 
gation , which conflicts with other studies (e.g jDe Wiin et ail 
120071) . While we cannot deduce a propagation velocity from 
the phase differences measured, we are confident that we in- 
deed detect u pward propaga t ion. 

Studies by ICenteno et al.l ([2006, 2009) show propagation 
in apparently vertical structures in a p l age re gion. Based 
on these observations, Khom enko et al.l d2008l) suggest that 
radiative relaxation is key to the propagation of 5-minute 
oscillations. Our analysis is based on obs ervations with 
much higher spatial resolution. Conceivably. ICenteno et alJ 
( 2006) observed slanted propagation inside their resolution el- 
ement, however, inversions of their spectropolarimetric mea- 
surements indicated that the field was truly vertical. In the 
results presented here no propagation is found in areas where 
the field can be reasonably expected to be vertical, suggesting 
strongly that inclination is critically important for the propa- 
gation of 5-minute oscillations, at least in this plage region. 

Given the results presented in this Letter, one may theorize 
that apparently vertical propagation is possible provided the 
field is sufficiently twisted. The field is then inclined with re- 
spect to the gravity vector, and oscillations may travel upward 
along this field. A structure with radius r reaches an inclina- 
tion of <f> at its periphery if it makes at least one complete twist 
over a height of h = 2-7rrtan(0). Since we do not accurately 
know the difference in the formation height of the observa- 
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tions used here, we cannot estimate the angle of propagation 
from these data. However, in order to adequately lower th e 
acoustic cut-off, 4> > 30° is needed dDe Pontieu et al.ll2004T) . 
Assuming a constant r = 50 km, we find h < 180 km. As the 
fluxtube expands with height, twist and hence inclination are 
enha nced naturally , so that h may be significantly increased 
(e.g.. lParkerlll974l) . Measuring twist in such a small structure 
is difficult, and so it is not immediately clear that this limiting 
value for h is reasonable. However, since we do not see signif- 
icant propagation inside the plage region, it seems that at least 
in this case, there is insufficient twist to allow for apparently 
vertical propagation of p-mode oscillations. 

4. CONCLUSIONS 

We have studied the propagation of waves in a small re- 
gion of plage using a Fourier phase-difference analysis. We 
find expected behavior of 3-minute oscillations. These oscil- 
lations are seen to propagate vertically, as is expected in the 
traditional view of the propagation of acoustic waves in the 
solar atmosphere. We find a much larger phase difference for 
these oscillations inside the plage region than outside it. Our 
results show that propagation of 5-minute p-mode oscillations 
happens only along inclined field at the periphery of the plage 
region. We find no propagation of p-modes in the core of the 
plage, where the field is expected to be mostly vertical. 
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